Skyrme pseudo-potential-based EDF parametrisation for 
spuriousity-free MR-EDF calculations 
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Abstract 

First exploratory steps towards a pseudo-potcntial-bascd Skyrme en- 
ergy density functional for spuriousity-free multi-reference calcula- 
tions are presented. A qualitatively acceptable fit can be accom- 
plished by adding simple three- and four-body contact terms to the 
standard central plus spin-orbit two-body terms. To achieve quanti- 
tative predictive power, higher-order terms, e.g. velocity-dependent 
three-body terms, will be required. 

1. Introduction 

Methods based on the use of effective energy density func- 
tionals (EDFs) [T] currently provide the only set of theo- 
retical tools that can be applied to all bound atomic nu- 
clei but the lightest ones in a systematic manner irrespec- 
tive of their mass, isospin, and deformation. One under- 
lying assumption is that of a universal effective energy 
functional that depends on normal and anomalous one- 
body density matrices and that re-sums the in-medium 
correlations whose details are irrelevant for low-energy nu- 
clear structure physics. Nuclear EDF methods coexist on 
two distinct levels. On the first level, that is tradition- 
ally called self-consistent mean-field theory, Hartree-Fock 
(HF) or Hartree-Fock-Bogoliubov (HFB), a single product 
state provides the normal and anomalous density matri- 
ces that enter the EDF. We call this type of method a 
single-reference (SR) approach. Exploiting the concept 
of symmetry breaking as a consequence of the nuclear 
Jahn- Teller effect [2], collective correlations that corre- 
spond to multipolc deformations or superfluidity can be 
easily modelled at the expense of losing good quantum 
numbers. On the second level, often called "beyond- mean- 
field methods", symmetry restoration and configuration 
mixing can be achieved within the Generator Coordinate 
Method (GCM) framework. Such method is referred to as 
a mult i- reference (MR) approach [3]. The MR techniques 
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aim at the explicit description of correlations that are re- 
lated to the small finite size of the nucleus and that neither 
can be easily absorbed into the EDF, nor be described by 
a single SR state. Typical examples are the mixing of dif- 
ferent shapes that coexist in a nucleus or the restoration 
of symmetries that are broken at the SR level. Besides de- 
scribing correlated ground states, the MR approach also 
gives access to excited states and transition moments be- 
tween them taking selection rules into account. 

Multi-rcfcrcncc techniques based on energy density 
functionals rely on the extension of the EDF to non- 
diagonal energy kernels. The generalised Wick theorem [T] 
provides the formal framework when the EDF is defined as 
the expectation value of a genuine Hamilton operator. The 
energy functionals widely used in nuclear physics, how- 
ever, are not of this form. We denote these as "general 
functional" in what follows. The standard procedure for 
the MR extension of such general functionals is made by 
formal analogy with the Hamiltonian case, cf. [5j |6] and 
references therein. It has been pointed out a few years ago 
that the usual procedure to set up the non-diagonal ker- 
nels of a general EDF is ill-defined |SJ H [SJ H ■ It gives 
rise to spurious contributions to the energy that can man- 
ifest themselves for example through divergences and/or 
finite steps when plotting the symmetry-restored energy as 
a function of a collective coordinate [H |S| ■ It also makes 
MR calculations return non-zero energies when restoring 
negative particle numbers [5J [TU] . These difficulties can be 
traced back to a breaking of the Pauli principle when set- 
ting up the EDF. In one way or the other, this is the case 
for all modern paramctrisations of the nuclear EDF. It is 
motivated either by phenomenology or for computational 
reasons [Tl ITT]. 

One popular choice for a general nuclear EDF is his- 
torically based on two- and three-body contact pseudo- 
potentials proposed by Skyrme [T2J. In the earliest ad- 
justments of its coupling constants for nuclear EDF cal- 
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culations it was used as a genuine Hamilton operatoiQ. 
This gave rise to three major problems in the resulting 
parametrisations: (i) the incompressibility of homo- 
geneous symmetric infinite nuclear matter (INM) was un- 
avoidably much larger than the empirical one; (ii) fits 
that constrain other INM properties such as isoscalar ef- 
fective mass mg/m, saturation density p sat and symme- 
try energy coefficient a sym to empirical values provided 
parametrisations that were unstable in the spin channel 
as indicated by the Landau parameter go being smaller 
than —1 |13j . and (iii) pairing matrix elements were far 
too small or even repulsive. Problem (ii) could be cir- 
cumvented reinterpreting the contact three-body force 
*3<5(ri — r 2 )<5(ri — r 3 ) as a density-dependent two-body 
force i h (l + x 3 Pf 2 ) Po I(n +r 2 )/2] <J(ri -r a ) with a = 1, 
where po(r) = p n (r) + p P (j) is the isoscalar density and 
P± 2 the spin-exchange operator, respectively, which al- 
ters the so-called "time-odd" part of the functional that 
governs the effective spin-spin interaction. Problem (i) 
could then be solved by reducing the exponent a from 
1 to values between about 1/6 and 1/3. And for most 
parametrisations of the Skyrme EDF, problem (iii) was 
removed by using different vertices in the particle-hole 
and particle-particle channels that were independently ad- 
justed to data. In addition, some specific terms in the 
EDF are often suppressed or approximated to simplify the 
numerical treatment [TJ. In the past this has led to high- 
precision parametrisations of the Skyrme EDF that kept 
a relatively simple form [TJ [14] and were numerically effi- 
cient, at the expense of sacrificing the exchange symmetry 
of the EDF. Other variants of the nuclear EDF have been 
set-up directly through their energy density without any 
reference to an underlying effective interaction (151 [US] . 

The breaking of the exchange symmetry that results 
from this practice introduces what is known as spuri- 
ous " self-interaction" in density functional theory for con- 
densed matter |17j . An analysis of the same problem in 
the context of effective interactions as used in nuclear 
physics has been given in Ref. [TB]. A similar spurious 
"self-pairing" appears when the normal and pairing parts 
of the EDF are not derived from the same effective interac- 
tion as pointed out in Ref. [8] ■ According to Ref. [B] , diver- 
gences and steps in MR EDF calculations originate from 
unphysical weights such self-interaction and self-pairing 
contributions arc multiplied with. Based on this analy- 
sis, a regularisation scheme that modifies those ill-defined 
weights while keeping the self-interaction itself has been 
proposed in Ref. [5] and applied with success to the case 
of pure particle-number projection in Ref. [5]. 

Not all general functionals are regularizable, though. 
The formalism can be applied only to functionals that 
correspond to polynomials in the density matrices [6j [9] . 
Parametrisations of a general regularizable Skyrme-type 
EDF of minimal form have been constructed recently [TT] 
and are currently used to test the regularised MR EDF 
scheme for general configuration mixing. Constraining the 
EDF to polynomial form makes its fit more difficult. 

In parallel, we also started the construction of function- 
als that are free of spuriosities from the outset by setting 

1 In the present context, so-called "density-dependent effective in- 
teractions" do not qualify as genuine Hamilton operators or "pseudo- 
potentials" . 



them up as the expectation value of a genuine Skyrme-type 
Hamilton operator taking all exchange and pairing terms 
into account such that the Pauli principle is obeyed. By 
construction, this removes all spurious contributions to 
the EDF at the price of having less independent terms 
in the EDF. In what follows, we denote such functionals 
a "pseudo-potential generated EDF". This replaces the 
problem of how to set up a formalism for the extension of 
a general and flexible EDF to non-diagonal kernels for MR 
calculations by the problems of what is the most efficient 
form for a predictive pscudo potential that can be straight- 
forwardly used in MR calculations and how to adjust it to 
data. This is a difficult task, as a pseudo-potential-bascd 
functional has much less independent coupling constants 
than a general functional of the same form. To be usable 
in SR and MR calculations incorporating pairing correla- 
tions, it is mandatory that the functional gives a reason- 
able description of the spin and pairing channels of the 
interaction. Otherwise the level sequence after angular- 
momentum projection is likely to be unrealistic. However, 
the difficulty to have even the correct sign of the inter- 
action in these two channels has been among of the ma- 
jor motivations to abandon Skyrme-type pscudo-potcntial- 
based EDFs in favour of the general ones in the 1970s. 

The aim of this contribution is to present the first ex- 
plorative steps towards the construction of a predictive 
pseudo-potential based EDF. We present the adjustment 
of a parametrisation that achieves an acceptable qualita- 
tive description in all channels of the interaction by adding 
three- and four-body contact terms without gradients to 
a standard two-body Skyrme operator. It will serve to 
benchmark MR EDF calculations and as a reference for 
what can be achieved without introducing derivatives in 
the new terms. 

2. Pseudo-potential-based EDF 

The pseudo-Hamiltonian used in the present work takes 
the form 

H = T« + «( 2Sk > + i/ 2C > + «( 3Sk > + z)( 4Sk > , (1) 

where J 1 *- 1 ' is the kinetic energy operator and v^ 2C ' the 
Coulomb interaction, which take their standard form. The 
{jl-WSk) are j\r_ DOC jy p ar ts of the Skyrme-type pseudo- 
potential for which we consider here the form 

£( 25fe ) = t (l + a; Pf 2 ) 5 Tl r 2 (2a) 

+ |(l + x iP£i) (k 12 2 V 2 + Kr^h) (2b) 
k 12 • 5 Ti r 2 k 12 (2c) 
+ iWo(o-i + d- 2 )-k^2 x 5 rir2 k 12 , (2d) 

V ^ ^ Vq ^ rir3 6 r2r ^ ^37*4 -\- S rir2 S r3r2 S r . 2r ^ -\- • • ■ ^ , (2f) 

where ti, Xi, Wo, uq and vq are unknown coupling con- 
stants, P±2 denotes the spin exchange operator, S r r is the 
Dirac distribution, <r is the vector of Pauli spin matrices 
and k 2 2 and k{| are the incoming and outgoing relative 
momenta. The three-body and four-body terms contain 
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three and twelve permutations of the coordinates, respec- 
tively. The SR and MR energy kernels are calculated as 



(PtShIp?? ,4**](r), (3) 



and nff 



where p** and re** denote normal and anomalous tran- 
sition density matrices, respectively. The resulting nuclear 
part of the energy density can be decomposed into bilinear, 
trilinear and quadrilinear parts according to their content 
in p and re 

£h = £ p tl 3 + £ft + £ P H PP + £ K H Kp 

+ £ P H PPP + £ h H Kpp + £% KKK . (4) 

The bilinear parts take the usual form [TH] 

£ B = E {^PlPl + + A* 1 * 1 ** ■ S q 

+ A S ^S q ■ S q + A TlPl TqPq + A TiP2 TqPg + A T ^T q ■ S q 

+ A^Tg ■ Sg + A^jg ■ j, + A^jg jg 

+ A V " lV ^ [V Pq ] ■ [V P g] + A V " lV ^ [V P g] ■ [V R g] 

I 4VS1VS2 [u 1 [\7 1 -U A J i J 2 J T_ 

+ A PlVJl P gV ■ Jg + A PlVJ2 P gV ■ Jg 



+ A^^ig ■ V X Sg + A^"j q ■ V X S A , (5 

fff = E l^' 1 ^^ + + A^fgP*g 

q ^ 



111' 



aJ*Ji J* J 1 /t^'^l 7* 7 

^2 ,J q.pp J q,vv ' ^3 ,J q,pv°q,vp 



(5b) 



They are functions of local matter density p q (r), kinetic 
matter density r 9 (r), current density j q (r), spin den- 
sity s 9 (r), spin- kinetic density T q (r), spin-current den- 
sity Jg lA1J /(r), pairing density p q (r), pairing- kinetic den- 
sity f 9 (r) and pairing-current density J q ,pv(j) |19j . Here, 
we choose a proton-neutron representation where the sum 
runs over neutron and proton densities q — n, p, and where 
q q denotes the other nucleon species. The trilinear 
parts of the energy density read 

£%» = E {A^pgpgp, + A™">B q ■ s qPq } , (6a) 
q 

£h P =Y. APlPlP2 ~PqhPq' (6b) 

q 

whereas the quadrilinear parts are given by 

£ PPPP = [AP^^PgPgPgPg + A^^Sg ■ S,S ? • Sg 

q 

+ AP^ S2 PgPgSg- ■ S^ , (7 &) 

£h PP = E { APlPlP2P2 P* q P q PqPq + A h ^p*~P q S q ■ Sg) , 

(7b) 



ckkkk _ AP1P1P2P2 X* p. n*h 

C H — l^i PqPqPqPq ■ 

q 



(7c) 



Each term in the functional is provided with a coefficient 
whose superscripts are labelled by numbers 1 and 2 that 
represent q and q, respectively. The coefficients are related 
to the parameters of the pseudo-potential @ through Ta- 
ble m 

3. A tentative fit 

^From the published work on pseudo-potcntial-based 
Skyrme EDFs from the 1970s it is clear from the out- 
set that the simple form of Eq. (flj will be unlikely to 
reach the quality of modern paramctrisations of the gen- 
eral Skyrme EDF. Our more modest aim is a paramctrisa- 
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Table 1: Coupling constants of the EDF expressed in 
terms of the parameters of the pseudo-potential. Missing 
entries are zero. 
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(MeV.fm 3 ) 
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Table 2: Coupling constants of SLyMRO. 

tion that could be used in time-reversal-breaking SR and 
angular- momentum restored MR EDF calculations based 
on HFB-type reference states. To accomplish this, we had 
to add several new constraints to our fit protocol. First, 
we had to avoid any unphysical instabilities, not just the 
ones related to too small values of the Landau parame- 
ters, but also finite-size instabilities of the kind discussed 
in Refs. [201 HI]- Second, we aimed at an overall repulsive 
interaction in the spin channels and at an attractive pair- 
ing interaction that provides gaps of realistic size. This, 
however, required to relax the constraints on many other 
properties, in particular those of nuclear matter. 

To construct an acceptable starting point for a fit, 
we took the parameters of SIV, which has a reasonable 
value of go and provides weakly attractive pairing, from 
Ref. [35] • We modified X\ and t\ to enhance pairing and 
then uq, vq and, to a lesser extend, £2 and x%, to bring the 
paramctrisation back to a better description of symmet- 
ric INM. Simultaneously, Xq had to be adapted to reject 
the onset of spin instabilities far enough above the satura- 
tion density. In a second step, the parameters were then 
fine-tuned to describe a set of nuclear masses and radii in 
addition to the nuclear matter properties in a least-square 
fit. There is no "best fit" in the usual sense, however. As 
several nuclear matter properties remain far from the em- 
pirical values, there are many possible fits of very similar 
(limited) quality. A typical parameter set obtained with 
that procedure that we will use for further studies is given 
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Table 3: Saturation density p sat , energy per particle E/A, 
symmetry energy a sym , effective mass raj, incomprcssibil- 
ity K^o and spin-Landau parameter go at saturation for 
the paramctrisations studied here. 
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Fig. 1: Neutron spectral gaps of singly magic even-even nuclei 
in the isotopic chains as indicated at the SR-EDF level. 



in Tab. [2] It is called SLyMRO in what follows. 

Table [3] compares some of the associated nuclear mat- 
ter properties with those of other paramctrisations used 
as pseudo-potentials in the literature. They differ in 
higher-order terms in Eq. ©. SV and SHZ2 arc pure 
two-body forces that have been used for isospin and 
angular-momentum projection of HF states (without pair- 
ing) [23] Hi]- SIII and SIV include the three-body term, 
whereas SLyMRO contains three- and four-body terms. 
Values obtained with a standard paramctrisation SLy4 of 
the general Skyrme EDF are given for comparison. 

As expected, all the pseudo-potential-based parametri- 
sations have difficulties to describe nuclear matter prop- 
erties. In particular, their effective mass is very low ex- 
cept for SIII. This leads to a density of single-particle lev- 
els around the Fermi energy that is much lower than the 
empirical one. For SIII, the larger effective mass leads, 
through the interrelations between the coupling constants, 
to spin instability as indicated by go < —1 [13] ■ For SIV, 
the reduction of the effective mass within the same func- 
tional form pushes go to values around zero. However, an 
analysis along the lines of Ref. [2Tj reveals that SIV has a 
finite-size spin instability in the S = 1, T = channel that 
cannot be revealed by the Landau parameters. Hence, nei- 
ther SIII nor SIV can be used in their pseudo-potential 
form in time-reversal invariancc breaking mean-field cal- 
culations. By contrast, such an analysis does not indicate 
unphysical instabilities for SV, SHZ2 and SLyMRO at den- 
sities relevant for low-energy nuclear structure. 

In Fig. [T] pairing properties are examined via the spec- 
tral gap -Epair, nj J d 3 r p n (r) of neutrons as obtained from 
spherical HFB-type SR calculations of singly magic even- 
even nuclei. Empirical pairing gaps obtained from a three- 
point mass difference are shown for comparison. When 
solving HFB equations, pairing matrix elements have been 
multiplied with a smooth cutoff at 8.5 MeV above and be- 
low the Fermi energy. SIII, SIV, SHZ2 and SV give null 
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Fig. 2: Binding energy residuals as a function of A for singly 
magic nuclei for spherical SR-EDF calculations. Nuclei in iso- 
tonic and isotopic chains are connected by lines. 



pairing or at least a weak pairing in some nuclei. Only 
SLyMRO, for which this property was enforced during the 
fit, provides pairing gaps of a realistic size. 

Figure [2] exhibits mass residuals for isotopic and isotonic 
chains of singly magic nuclei as obtained from spherical 
HFB-typc SR calculations. The particularly large drift of 
the curves obtained with SLyMRO results mainly from its 
very low value for the asymmetry energy coefficient a sym , 
which cannot be increased without jeopardising pairing or 
the time-odd terms in the EDF. We have checked that, in 
spite of its poor description of masses, SLyMRO gives a 
reasonable description of the deformation of nuclei in the 
sd and pf shell regions, and of their rotational bands in 
cranked HFB calculations. 

4. Conclusion 

The simple pseudo-potential ^ allows for a description 
of basic nuclear properties that does not meet the stan- 
dard of state-of-the-art SR-EDF calculations. However, it 
is sufficient for its main purpose of benchmarking multi- 
dimensional MR EDF calculations, which then will be free 
from the pathologies encountered with any modern stan- 
dard parametrisation. In particular, we achieved a fair 
description of pair correlations while avoiding (finite-size) 
spin-instabilities. However, the simple form used here is 
clearly over-constrained and higher-order terms will be 
needed to replicate the performance of the best available 
general Skyrmc EDFs at the SR level. Work in that direc- 
tion is underway. As a first step, the most general central 
three-body contact operator up to second order in gradi- 
ents has been worked out recently [33] . First fits that aim 
at the set-up of a suitable protocol are currently underway 
and show promising results. 
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